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The small-world transition is a first-order transition at 
zero density p of shortcuts, whereby the normalized shortest- 
path distance £ —J/L undergoes a discontinuity in the ther- 
modynamic limit. On finite systems the apparent transition 
is shifted by Ap ~ L~ d . Equivalently a "persistence size" 
L* ~ p _1//<i can be defined in connection with finite-size ef- 
fects. Assuming L* ~ p T , simple rescaling arguments imply 
that r = 1/d. We confirm this result by extensive numerical 
simulation in one to four dimensions, and argue that r = 1/d 
implies that this transition is first-order. 

A physically meaningful characterisation of disordered 
connection networks is given by their average chemi- 
cal distance or shortest-path distance £ = l/NJ^.-y < 
£ij >, where £ t j is the minimum number of links that 
must be traversed to join sites i and j, and <> means av- 
erage over disorder realizations. If the network connectiv- 
ity is topologically d-dimensional, I is proportional to L, 
its linear dimension. But on randomly bonded networks 
containing L d sites one has £ ~ log(L) much shorter than 
L. Small- world networks [Q interpolate between these 
two cases. They consist of a regular d-dimensional lattice 
with a small fraction p of long-range bonds, or shortcuts, 
per site (i.e. there is a total of B = pL d long-range bonds 
on top of a regular lattice). These systems have received 
attention recently , following the observation jjj that 
a very small fraction of long-range bonds is enough to 
produce a transition to a regime where £ ~ log(L). 

Barthelemy and Amaral (in the following BA) [Q re- 
cently showed that for any nonzero p, the logarithmic 
regime is always reached if large enough systems are con- 
sidered. In other words there is a "crossover size" L*(p) 
such that for system size L < L* , £(L,p) scales linearly 
with L, while for L > L* £{L,p) grows as logL. This 
behaviour is described by the scaling form 

l(L,p) ~ Lf{L/L*(jp)) (1) 

with / a scaling function such that fix) — > constant for 
x « 1 and / — > log(x)/x for x » 1. BA however claim 
that this "small world transition" is not a phase transi- 
tion but a "crossover phenomenon" . These authors fur- 
ther propose that L*(p) ~ p~ T for p — > 0. By numerical 
analysis on relatively small systems, they find r w 2/3 in 
one dimension and conclude that r < 1 provides evidence 
for the importance of "nonlinear effects" . 

Some of their conclusions have been subject to crit- 
icism recently. Firstly Barrat jq] stressed the obvious 
fact that r cannot be smaller than 1 in one dimension, 
otherwise the total number of long-range bonds at the 
"crossover point" would go to zero with system size. Bar- 
rat confirmed by numerical simulation on large systems 



that T\d — 1. Similar arguments can be used to conclude 
that r > 1/d in d dimensions. 

A second controversial point concerns the nature of 
the small-world transition, which is the subject of this 
letter. Newman and Watts (thereafter NW) recently 
suggested that the transition from logarithmic to linear 
scaling of £ is a normal second-order phase transition 
at p — 0, i.e. one with a diverging characteristic length 
~ P ~ v i an d derive v = 1/d by renormalisation-group 
(RG) analysis. 

In this Letter we show that the small-world transition 
is a first-order transition at p = 0, where the "order pa- 
rameter" C(p) — £/L undergoes a discontinuity in the 
thermodynamic limit. Within this picture, L* ~ p~ T in 
([!]) is a persistence size associated with finite-size correc- 
tions. We derive r = 1/d by rescaling arguments similar 
to those of NW and confirm this result by extensive nu- 
merical simulation in one to four dimensions. The fact 
that L* ~ p _1 / d implies that, on finite systems, there is 
a shift or broadening of order Ap ~ L~ d in the apparent 
transition point. In support of our claim for a first-order 
phase transition we recall earlier work of Nicnhuis and 
Nauenberg H and of Fisher and Berker JhJ, who first 
discussed the relationship between first-order transitions 
and finite-size corrections of order L~ d . 

Let us in the first place show that, if Eq. ([!]) holds with 
L* ~ p~ T , then r = 1/d. The first steps of our deriva- 
tion are similar to those in NW ||. Consider a decima- 
tion transformation with rescaling parameter b << L, 
i.e. define "blocks" of linear size b, each containing b d 
sites of the original lattice. After rescaling we are left 
with a lattice of linear dimension L = L/b. Since all 
linear dimensions are rescaled by b, one also has that 
L* = L* /b. Two blocks I and J in the decimated lattice 
are connected by a long-range bond if and only if any 
pair of original sites i <Z I and j C J are connected by a 
long-range bond. Since the total number B of long-range 
bonds is invariant under rescaling (because their typical 
length is ~ L, much larger than b), one has pL d = pL d , 
and therefore p — pb d . After rescaling one must have 
L* — p~ T and therefore b = L*/L* = (p/p)~ T = b dT , 
which implies t = 1/d. 

Our extensive numerical investigations confirm this 
theoretical prediction. We used hypercubic lattices with 
helicoidal boundary conditions, of sizes up to 250000 lat- 
tice sites in one to four dimensions, to which random 
bonds were added with probability p per site. Since 
in this problem all bond costs are unity, shortest-path 
distances £ji can be determined efficiently by means of 
Breadth-First-Search (BFS) whose time-complexity 
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is O(N), and is simple to program. We calculate £(L,p) 
by averaging over 1000 — 2000 disorder realizations for 
several lattice sizes and p values. From Eq. ([!]), the 
asymptotic derivative dl/dlogL becomes proportional 
to L* as L — ► oo. We thus fit a straight line to the J 
vs. \ogL curve, using the four largest sizes simulated for 
each value of p, and obtain the values of L* (p) displayed 
in Fig. [j]. Fitting these data we obtain r = 1/d within 
5% error in all cases. 
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FIG. 1. Persistence size L* vs p in one to four dimensions. 
Our fits for r such that L* ~ p~ T agree, within 5% error 
in all cases, with r = 1/d, the value predicted by reseating 
arguments. 

Now consider the order parameter 
C(p) = limL^oo^/i, which is zero for all p > 0, but 
a (lattice-dependent) constant for p = 0. This suggests 
that it is reasonable to identify the small-world transi- 
tion as a discontinuous, or first-order phase-transition at 
p = 0. Our characterisation of the transition as first- 
order is however not based on the discontinuity of C, but 
a consequence of having t = 1/d as we now discuss. 

As first shown by Nienhuis and Nauenberg M in the 
context of the Renormalisation Group, and later more 
generally by Fisher and Berker (To) , at a first-order fixed 
point one eigen-exponent must take the value y = d 
(this is necessary in order to have phase-coexistence) and 
therefore finite-size corrections of order L~ d can be ex- 
pected. Conversely, if y = d then some first-derivative 
of the free-energy density undergoes a discontinuity at 
the transition, i.e. y = d implies that the transition is 
first-order. The existence of an eigen-exponent y = d 
gives rise to finite-size corrections (e.g. the shift in the 
apparent transition point) which are of order L~ d p0[ . 
We see then that finite-size corrections of order L~ d are 
a signature of first-order transitions. 

For pedagogical reasons let us first review how finite- 



size corrections of order L~ d arise at a typical first- 
order transition, and later make the analogy with small- 
world transitions explicit. Consider an Ising-type d- 
dimensional ferromagnet of linear size L at very low tem- 
perature, subject to an external magnetic field h. Assume 
that a single spin (e.g. the central one) is "pinned down" . 
This system undergoes a first-order phase transition at 
h = 0. In the thermodynamic limit, the magnetisation 
is m(h) « +1 for h > 0, and m(h) « —1 for h < 0. 
For L finite and h > 0, most of the spins will point up 
if the bulk energy hL d associated with the field is larger 
than the energy Jj associated with breaking the 7 bonds 
around the pinned spin. Therefore on finite systems of 
size L, the transition no longer happens at h = but at 
hL d — Jj = 0, i.e. it suffers a shift Ah of order L~ d . 
Equivalently, for each finite h — > + a persistence size 
L* ~ hr x l d can be defined such that m < for L < L* , 
while m > for L > L*. 

In the case of small- world systems, the order parameter 
C{p) is a step function at p = when L — > 00. But for 
finite sizes L, C(p,L) — £(p,L)/L is nonzero not only 
at p = but for all p < Ap ~ L~ d , or equivalently for 
L < L* = p~ x / d when p is fixed. The situation is clearly 
the same as in our previous example, i.e. the transition 
suffers a shift of order L~ d on finite systems. 

Thus L* is not a "crossover length" S but reflects 
the existence of boundary-dependent finite-size correc- 
tions around a first-order transition. We thus prefer to 
call L* "persistence size" JlC| ] . 

As discussed in the previous paragraphs, by showing 
that the finite-size corrections exponent r equals 1 / d we 
have also demonstrated that the small- world transition is 
first-order. But contrary to what is sometimes assumed, 
the first-order character of a transition does not preclude 
the existence of critical behaviour. Notice that C behaves 
in the same way as the magnetisation density m(T) of 
a one-dimensional Ising model, or the spanning-cluster 
density P^i^p) in a one-dimensional percolation model. 
These two cases do display critical behaviour, and consti- 
tute simple examples of first- order critical points, respec- 
tively at T — and p = 1. First-order critical points have 
been briefly discussed by Fisher and Berker (FB) | fL0f . 

It is usual (though not entirely correct) to refer to Id 
percolation or Id Ising as displaying "second-order" tran- 
sitions. The correct designation is "first-order critical". 
Around a FOCP finite size corrections Ap are of order 
L~ d because this is true for any discontinuous transition. 
Since at a critical point one also has Ap ~ L~ x l v this in 
turn means that v = 1/d for a FOCP pQ] . Thus the fact 
that the global observable £ for finite small-world net- 
works is a function of a single scaling variable x = Lp x / d 
is compatible both with a non-critical first-order transi- 
tion (in which case x = L/L* with L* a persistence size) 
and with a first-order critical point (in this case x — L/£ 
with £ ~ p^ 1 ^ an internal characteristic length). 

Newman and Watts HJfl characterise the small-world 
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transition as a "normal second-order" phase transition, 
i.e. they assume critical behaviour ("first-order critical" 
would be the appropriate term if there were critical be- 
haviour). As said, the finite size behaviour of £(L,p) 
cannot be used to distinguish between a normal first- 
order point and a FOCP. A possible way to prove or 
disprove the existence of critical behaviour in small-world 
networks is to look at bilocal observables, and see whether 
an internal correlation length can be identified that gov- 
erns the physics of the transition and diverges at p c . We 
study in this work the behaviour of a particular bilocal 
observable: £(r), the average shortest-path distance be- 
tween two points separated by an Euclidean distance r. 
The basic question we wish to answer is whether the be- 
haviour of £(r) in the thermodynamic limit is dictated by 
a finite characteristic length £(p) that diverges at p = 0. 
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FIG. 2. We divide the system into blocks (shaded) 
of size r (the Euclidean distance between a and b) in or- 
der to estimate the average shortest-path distance £(r) when 
r is much larger than the mean separation p~ 1//d between 
shortcut-ends (white dots). By considering the distance trav- 
elled on the lattice (dashed lines) along any path from a to b 
through shortcuts (black lines), one finds that £(r) — r unless 
r > fT 1/d log L (see text). 

Newman and Watts Q recently argued that the char- 
acteristic length £ for the small-world transition is the 
mean separation s(p) ~ p~ x l d between shortcut-ends. 
Their interpretation is consistent but they apparently 
base this claim on the sole fact that the global observ- 
able £ is a function of L/s(p) only. No analysis is done 
of any bilocal observable, which we think is the only way 
to assess the existence of critical behaviour. 

We next show that s(p) is not a lengthscale dictating 
the behaviour of shortest-path lengths £(r). We will find 
that a characteristic length r c does in fact exist for £(r), 
but that this characteristic length is r c ~ s(p) logi, i.e. 
it diverges with system size, for all p > 0. 

Consider the shortest-path distance £{r) between two 



points a and b separated by an Euclidean distance r. If 
r << s(p) = p~ x l d then clearly £(r) = r. Let us then 
discuss the case s(p) << r << L. Imagine dividing the 
lattice into blocks of linear size r as shown in Fig. ||. 
The resulting rescaled system can be viewed as a ran- 
dom graph made up of N = (L/r) d nodes, randomly 
connected by a total of B = pL d long-range bonds. The 
average coordination number of a node (the number of 
shortcuts connected to it) is k ~ pr d . 

Consider now a path from a to b through one or more 
shortcuts (Fig. |^). In order to determine the "cost" of 
such path, we will associate a cost zero with each short- 
cut, and only count the total distance travelled on the 
Euclidean lattice within each visited node (dashed lines 
in Fig. U). There are typically no shortcuts with both 
ends in the same block because r << L. Thus a path 
(a, b) through shortcuts must contain two or more out- 
going shortcuts and visit one or more extra nodes. Thus 
we are looking for closed "loops" in the random graph, 
that start and end in the original node (the one that 
contains a and b). Now we know that, starting from any 
given node, after n shortcut-steps a total of ~ k n nodes 
can be visited in average. A closed loop will exist if one 
of these visited nodes is the original one. The probabil- 
ity for this to happen is simply k n /N, which is negligi- 
ble unless n ~ logiV/ \ogpr d . Thus the shortest closed 
loop in the random graph of blocks will typically contain 
O(log N/ \ogpr d ) or more nodes. The cost associated 
with visiting a node can be bounded to be at most 0(r) 
and at least 0(s(p)). Taking a conservative estimate let 
us say it is s(p). Then the lowest cost of going from a 
to b through shortcuts is typically ~ s(p) log N/ log pr d . 
Thus it will be "cheaper" to travel from a to b directly on 
the Euclidean lattice than through shortcuts (and conse- 
quently £{r) = r) if r < s(p) log Nr~ d / log pr d . 

This approximate scaling argument is confirmed in its 
main aspects by recent analytical Q and numerical 
studies (see Fig. ||), which show that £{r) behaves as 
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r ioi- r <r c =p- 1/d log{p 1/d L) 
r c for r > r r 



(2) 



In other words, the characteristic length for shortest- 
path distances £(r) is not s(p) but r c above, which di- 
verges with L for all p. 
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a non- critical first order transition despite the fact that 
s(p) diverges at p = 0. We feel that a similar phenomenon 
happens in the small-world transition since the average 
distance between shortcut ends s(p) diverges at p = 
but plays no role in determining shortest paths. 
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FIG. 3. Shortest path distance i(r) rescaled by r c 

vs. the rescaled Euclidean distance r/r c , as measured on 
one- dimensional small-world lattices [[jj. The characteristic 
length r c ~ p~ 1 ^ d \ogL diverges with system size. 

In summary, we showed that the small- world transition 
is neither a "crossover phenomenon" Q nor a second- 
order phase transition || , but a first-order phase transi- 
tion at p = 0. Finite size corrections are dictated by a 
"persistence size" L*(p) which behaves as p~ T ||. Sim- 
ple rescaling arguments similar to those used by Newman 
and Watts |J show that r = l/d. This is confirmed by 
extensive numerical simulations in one to four dimen- 
sions. We argue that having r = l/d implies that the 
transition at p = is first-order P,|i"o|], and discuss the 
possibility of further characterising this first-order point 
as a first- order critical point pofl . This can only be done 
by looking at bilocal observables and we choose to study 
£(r). A relatively simple scaling argument, confirmed 
by more rigorous studies shows that the character- 
istic length r c dictating the behaviour of £(r) diverges 
with system size, i.e. it is not possible to define a fi- 
nite p-dependent characteristic length for £(r) in the 
thermodynamic limit. We suggest that s(p) as defined 
by NW U is an "irrelevant" internal length, i.e. one 
diverging at p = but playing no role in determining 
shortest paths. A simple example of such an irrelevant 
length is given by the following example: consider fcr- 
romagnet at T << T Cl h = and with one site "pinned 
down" . Now instead of a magnetic field h, assume there 
is a density p of randomly located sites whose spins are 
"pinned-up" . As in our previous example, a first-order 
transition happens at p = where the magnetization m 
changes sign, and this transition is clearly non-critical. 
In this "tailored" example, an internal characteristic size 
s(p) ~ p- x l d (the average distance between pinning sites) 
can be defined, which however plays no role in the physics 
of the magnetic transition. This magnetic transition is 
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